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We study the charge ordering (CO) in the one-dimensional (ID) extended Hubbard model at 
quarter filling where the nearest-neighbor Coulomb repulsion and dimerization in the hopping pa- 
rameters are included. Using the cluster mean-field approximation to take into account the effect of 
quantum fluctuations, we determine the CO phase boundary of the model in the parameter space at 
T = K. We thus find that the dimerization suppresses the stability of the CO phase strongly, and in 
consequence, the realistic parameter values for quasi-lD organic materials such as (TMTTF^PFe 
are outside the region of CO. We suggest that the long-range Coulomb interaction between the 
chains should persist to stabilize the CO phase. 
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I. INTRODUCTION 

The charge-ordering (CO) phase transition has re- 
cently attracted much attention in physics of strongly 
correlated electron systems such as transition-metal ox- 
ides and organic charge-transfer salts. A possible mech- 
anism of CO is the localization of electrons on a lattice 
due to long-range Coulomb repulsions, which is a lattice 
version of the Wigner crystallization. The spin degrees 
of freedom of the system of CO are yet active and can 
undergo an additional phase transition by lowering tem- 
perature, to result in a variety of magnetic ground states 
such as spin-Peierls (SP), spin-density-wave (SDW), and 
antiferromagnetic (AF) states. One of the simplest mod- 
els that allow for such CO is the one-dimensional (ID) 
Hubbard model with onsite (U) and nearest-neighbor 
(V) Coulomb repulsions, where the charge density wave 
(CDW) with two-fold periodicity, i.e., 4fcp-CDW with kp 
being the Fermi momentum, is realized when the band 
filling is either 3/4 or 1/4 (quarter-filling) . 

Bechgaard salts (TMTSF^X and their sulfur analogs 
(TMTTF) 2 X, where X is PF 6 , Br, C10 4 , etc., offer a se- 
ries of materials suitable for studying quasi- ID correlated 
electron systems at quarter filling p). It is known that 
the lattice dimerization of the systems causes an alterna- 
tion of hopping integrals and makes the system a Mott 
insulator. The metal-insulator transition of the series is 
then controlled by the interchain hopping parameter t±, 
i.e., by the dimensional crossover from ID to 2D j|, ||, Q. 

Recently, clear evidence of the CO phase transition in 
(TMTTF^PFg has been given by the measurements of 
dielectric response || and NMR spectroscopy Q. This 
material has the smallest interchain coupling t± in the se- 
ries and is in the ID confinement regime 0, §. The tran- 



sition temperature of CO is reported to be Tco — 100 
K, which is very high for small energy scales of organic 
systems. A new research area for studying CO has thus 
been established 

A number of theoretical studies have so far been made 
on the ID extended Hubbard model at quarter filling 
@, |, § [TTJ, [l| If, 0, ||? ], which provide useful 
information to consider the physics of CO. While some of 
them have not included the effect of lattice dimerization 
1^0^, recent Id calculations p3|, 
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have taken into account the effect and explained 
how the 2kp- or 4fcp-CDW coexists with 2fcp-SDW in 
the ground state of the systems. We note however that 
the mean-field analyses usually overestimate the stability 
of the ordered phase. Coupling to the lattice degrees of 
freedom has also been argued to play an essential role 

0- 

Motivated by such development in the field, we study 
in this paper the CO of the ID extended Hubbard model 
with the dimerization of hopping parameters, and con- 
sider its implication to the observed CO in the TMTTF 
family of organic conductors. We use a new method 
of calculating the long-range order in strongly-correlated 
electron models, i.e., an extension of the so-called clus- 
ter mean-field approach jl8| to correlated fermion sys- 
tems. Thus, the amplitude of CO, rather than its corre- 
lation functions or exponents, can be calculated directly 
as in the standard mean-field theory. Moreover, criti- 
cal interaction strength of the quantum phase transition 
can be evaluated with improved accuracy since the ef- 
fect of quantum fluctuations is automatically taken into 
account. The CO phase boundary and its dimerization 
dependence of the model is thus determined in the pa- 
rameter space (U, V) at T = K. 

We will thereby show that the critical interaction 
strength V c is much larger than the value obtained in 
the mean-field approximation, as expected, and also that 
the dimerization in the hopping integrals suppresses the 
stability of the CO phase strongly. Then, it follows that 
the present model does not have the CO ground state 
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FIG. 1: Schematic representations of (a) the Hubbard chain 
and (b) the chain with the mean-field bond. 

if we assume realistic values of the electronic parame- 
ters for quasi-lD organic materials (TMTTF) 2 PF 6 and 
(TMTSF) 2 C10 4 . This result suggests that the present 
framework of the model under the influence of weak 
three-dimensionality (3D) is not sufficient for describing 
the CO observed in (TMTTF) 2 PF 6 . Inclusion of any 
additional degrees of freedom, in particular the 3D long- 
range Coulomb interaction between the chains, may be 
essential for stabilizing the CO phase. 

II. MODEL AND METHOD 

The extended Hubbard model is defined by the Hamil- 
tonian 

L/2 

H= - ti£ $>Li,.C2vr + H.c) 

i=l a 
L/2-1 

i—1 a 

L L-1 

+ U^2n ii - t n ii i + V'^2n i n i+1 (1) 

i=l i=l 

on the ID lattice of size L (even), where c\ a (cj iCT ) is the 
electron creation (annihilation) operator at site i and spin 
a (=|, |), and rii = n^f + is the number operator. 
We introduce a dimerization in the hopping parameters, 
ti > i 2 . U and V are the onsite and nearest-neighbor 
Coulomb repulsions, respectively. We restrict ourselves 
to the case of quarter-filling. 

Our cluster mean-field method is the following: We 
first introduce an exactly solved finite-size system (or a 
cluster, see Fig. 1) and then assume the presence of the 
mean-field acting on the edges of the system. The value 
of the mean-field is then determined self-consistently so 
as to minimize the total energy of the system. In the 
present case, we prepare an L-site chain treated exactly 
by a numerical method and assume that the mean-field is 
applied on the edge sites 1 and L of the chain. We then 
replace the mean-field by the 'mean-field bond' connect- 
ing between the sites 1 and L; i.e., we assume 

nxriL ~ ni(n L ) + {n\)n L - (ni)(n L ) (2) 
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FIG. 2: Converged charge-density distributions for the models 
without dimerization (t\=t2=t). The cluster of length L = 8 
is used. In (a), three typical solutions are shown: type I 
(squares), type II (circles), and type III (triangles), which are 
discussed in the main text. In (b), oscillations of the charge 
density near the critical point at U/t = 9 are shown. 



for the bond, where (ni) and {njy) are the mean-fields at 
sites i = 1 and L. The Hamiltonian of the bond 

H\L = -t 2 y^Ci.gCL.o- + H.c) 
a 

+ V{ ni {n L ) + (m)n L - {m)(n L )), (3) 

is then added to the Hamiltonian of the cluster Eq. (1). 
The total Hamiltonian is diagonalized numerically by the 
Lanczos technique on small clusters, and the mean-fields 
are evaluated as the expectation values of n\ and til for 
the ground state. The iterations are made to achieve 
self-consistency in the values of the mean- fields {n\) and 
(ul)- Converged solutions are obtained after 30 — 70 it- 
erations depending on various conditions including the 
initial values of (ni) and (til)- We need to try a number 
of the initial values in order to confirm that the iteration 
converges to the unique lowest-energy solution. In some 
cases, we find two (or possibly more) different solutions, 
of which we choose the lowest-energy one; we should note 
that the level crossing in the lowest-energy states some- 
times occurs in small-size systems (which is apparent, 
e.g., in Fig. 3(a) given below). 

In Fig. 2(a), we show the converged charge distribu- 
tions for three typical cases: the type I is for the paramag- 
netic metallic state where a uniform charge distribution 
is seen, the type II is for the CO state where the oscilla- 
tion of two-fold periodicity is clearly seen, and the type 
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III is for the par ameter region corresponding to phase 
separation Jj| [Ll|. In Fig. 2(b), we show how the CO 
oscillation ceases near the critical point with decreasing 
V; we find that the 4kp oscillation becomes ill-defined 
between the interaction strength V/t = 2.67 and 2.68, 
at which we decide the critical interaction strength V c is 
located. These are the cases without dimerization but 
the situation is similar even when the dimerization is in- 
troduced. We should note that, in some cases, a 2/cf 
oscillation (corresponding to the effect of Friedel impu- 
rity scattering) is superimposed on the 4/cp CO oscilla- 
tion, which is however rather small except in the region 
of phase separation. 



III. RESULTS AND DISCUSSION 

The obtained CO phase boundary of the model is 
shown in Fig. 3. In case where there is no dimeriza- 
tion (ii = t% = t), we can compare our results with 
those of the previous studies [0, g], where the charge 
gap and Luttinger- liquid parameter K p were evaluated in 
small-size systems and the metal-insulator phase bound- 
ary (which coincides with the CO phase boundary when 
ti = £2) was determined. We find that our results are 
in good agreement with the previous results. Exact so- 
lutions are available for two limiting cases 0: V c /t = 2 
at U — » 00 and U c /t — 4 at V — > 00, which are also con- 
sistent with our results. We note that a rather inhomo- 
geneous charge distribution is obtained in the parameter 
region V S> U as seen in Fig. 2(a); in this region, the 
model is reported to have exotic phases where the su- 
perconducting pairing fluctuations are dominant or the 
phase separation occurs |)| [Tl| . 

The cluster-size L dependence of the critical interac- 
tion strength V c for t\ = £2 is shown in Fig. 4. An os- 
cillatory behavior corresponding to the so-called shell ef- 
fect for periodic boundary condition is found between 
L = 8m (open shell) and L = 8m — 4 (closed shell) 
where m = 1,2, • • •. We find that the values of V c thus 
calculated for L — 8, 12, and 16 appear to be extrapo- 
lated well to the value at 1/L — > which is estimated 
from the density-matrix renormalization-group (DMRG) 
calculation of the charge gap A c : The results of A c at 
1/L — > for several values of V are fitted to the expo- 
nential dependence jl^| A c /£ — aexp[— (3t/(V — V c )] (ob- 
tained from weak-coupling theory near the critical point), 
from which we determine the value of V c (in the accuracy 
of ±0.05£). The details of our DMRG calculation is given 
in Ref.|£0| We expect from these results that the value 
of V c at t\ > £2 also should not have too strong size- 
dependence when 1/L — ► 0. 

In case where the dimerization is present (t\ > £2), 
we first note that the charge gap opens Q except for 
the noninteracting case (U = 0); an infinitesimal U value 
is enough to make the Umklapp process relevant in the 
renormalization, leading to the opening of the charge gap. 
When U and V becomes large, we find that the quantum 
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FIG. 3: Dimerization ti/t2 dependence of the CO phase 
boundary. Results are shown for t\/t% = 1 (squares), t\/ti = 
1.11 (circles), and t\/ti = 1.43 (triangles). Clusters of the size 
(a) L — 8, (b) L = 12, and (c) L = 16 are used. Upper-right 
side of the boundary is the CO phase, and lower-left side of 
the boundary is either paramagnetic metallic phase (t\ = £2) 
or Mott insulating phase (ti > £2). 

phase transition from this Mott insulating state to a CO 
state occurs. We can detect this transition by observing 
the appearance of oscillations in the calculated charge 
distribution. Note that the mean-field bond is assumed 
to have the hopping parameter ti which is smaller than 
£1, i.e., the bond is chosen so as to connect two 'dimers' 
(where the dimer is a pair of sites with the larger hopping 
parameter). This is because, in the strong dimerization 
limit, only a single electron is present in a dimer, and 
thus V should not work in a dimer. 

In Fig. 3, we find that the CO phase boundary ex- 
hibits remarkably strong dimerization dependence; i.e., 
the dimerization suppresses the stability of the CO phase. 
General features are the following: (i) with increasing 
dimerization, the boundary in the large U region shifts 
to the left, i.e., V c /t2 becomes large, and (ii) the bound- 
ary in the large V region shows a small upward shift 
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FIG. 4: Cluster-size L dependence of the critical interaction 
strength V c . Results are shown at U/t = 7, 9, and 12 for the 
model without dimerization (ti = ti = i) . The values of V c 
at 1/L — > determined from the calculation of the charge gap 
by the density-matrix renormalization group method are also 
shown for comparison. 



with increasing dimerization. These features can be un- 
derstood as follows: We first note that, in the limit of 
strong dimerization t\ » t2, the critical V value is de- 
termined as V c = At\ when U — > oo. This is because 
in this limit the charge (which is either on the left or 
on the right site of a dimer) may be expressed as the 
pseudospin and the system becomes equivalent to the 
ID quantum Ising model p2f , the critical point of which 
is determined by the competition between the quantum 
fluctuation t\ of an electron in a dimer and the interac- 
tion V between the two dimers leading to CO. We there- 
fore have 14/(4^2) = t\/t2 — > 00 when - * 00, which 
explains the behavior (i). When V — > 00 and U is fi- 
nite, we find from an estimation of kinetic energy of the 
lowest-energy charge excitation that the critical U value 
is given as U c = 2(<i + t 2 ), irrespective of the strength of 
dimerization. Thus, we have U c /t2 = 2{t\/t2 + 1), which 
gives a contribution to the behavior (ii). A clear upward 
shift of the phase boundary is however not seen in Fig. 3, 
which is partly because sufficiently large values of V and 
t\ jt2 are not used and partly because the finite-size effect 
becomes apparent here. 

Besides the limiting cases discussed above, we can com- 
pare our results with the results of the mean-field cal- 
culation pUfl . Our results shown in Fig. 3 are for two 
dimerization strengths, £2/^1 = 0.9 and 0.7, which cor- 
respond to the realistic values for (TMTSF)2C104 and 
(TMTTF) 2 PF 6 , respectively. We then point out that the 
critical interaction strength obtained here is much larger 
than that of the mean-field calculation: e.g., we obtain 
V c /t 2 ~ 3 at U/t2 = 5 and £2/^1 = 0.9, while the mean- 
field calculation jl3| gives V c /t 2 — 0.65. This is also the 
case without dimerization (t\ = £2 = t): e.g., we obtain 
V c /t = 2.9 at U/t = 5, while the mean-field calculation 
|fl3| gives V c /t = 0.4. These discrepancies stem from the 
strong quantum fluctuations of the present model. 

Finally, let us consider experimental implication of our 



results. A recent careful estimation [£0| of the values of 
the parameter V for real materials has given V = 0.21 
eV (or V/t 2 = 0.8) for (TMTSF) 2 C10 4 and V = 0.18 
eV (or V/t 2 = 2.0) for (TMTTF) 2 PF 6 . Also the val- 
ues U/t 2 = 5.6 for (TMTSF) 2 C10 4 and U/t 2 = 10 for 
(TMTTF) 2 PF 6 have been reported @, These val- 
ues of the parameters in the (U, V) plane (see Fig. 3) are 
thus located in the Mott-insulating uniform phase, far 
apart from the CO phase boundary. We therefore con- 
clude that the strength of the nearest-neighbor Coulomb 
repulsion for these two materials is too small for the cor- 
responding ID extended Hubbard model with dimeriza- 
tion to be in the region of CO. Since there is no CO at 
T = K, it seems quite unlikely that the present model at 
T > can have the CO even on the implicit assumption 
of the presence of weak 3D coupling. Then, to explain 
the CO observed in (TMTTF^PFg, it seems necessary 
to take into account any additional degrees of freedom 
in our ID model; a comparatively large value of Tqo 
may suggest that there exists a rather strong interchain 
coupling via the long-range Coulomb interaction, which 
works to stabilize the CO phase. We may say that, while 
the electron conduction of (TMTTF) 2 PF 6 at T > T co 
is of the ID nature |Q, 0], the observed CO is due to the 
presence of the 2D-3D Coulombic coupling between the 
chains. Coupling to the anionic potential may also play 
an important role in the stabilization of CO as has re- 
cently been suggested in Ref. (23). Whether these effects 
may be renormalized into the parameter values of our ID 
model is a highly nontrivial problem; the issue is thus left 
for future work. 



IV. CONCLUSION 

We have studied the CO in the ID extended Hubbard 
model at quarter filling by introducing the dimerization 
in the hopping parameters as well as the nearest-neighbor 
Coulomb repulsion. We have used the cluster mean-field 
approximation to take into account the effect of quan- 
tum fluctuations, and have determined the CO phase 
boundary and its dimerization dependence of the model 
in the (U, V) plane at T = K. We have thus found that 
the realistic parameter values for (TMTTF^PFg and 
(TMTSF) 2 C10 4 are well outside the region of CO. We 
have argued that the present ID model does not provide 
a sufficient framework for describing the CO observed in 
(TMTTF) 2 PF 6 . 
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